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Abstract. I review recent results from lattice QCD calculations on the in-medium hadron 
properties. I discuss the thermal dilepton rates, heavy quarkonium properties as well as 
the chiral and 1/(1)a symmetries at finite temperature. 



1 Introduction 

Understanding the phases and properties of strongly interacting matter at high temperature and density 
has been one of the central goals of research in nuclear science. The hadrons inside the medium would 
feel the modifications when the medium is heated or compressed [ 1 ]. The vector meson p and the axial 
vector meson a\ are degenerate when the chiral symmetry is restored from the spontaneously broken 
phase at a certain temperature. The particles 77' and 77 are degenerate when the U(1)a symmetry is 
restored. The heavy quarkonia, e.g. J/tf/ and T, due to their small size and large mass, may survive up to 
certain critical temperatures and could serve as medium thermometers 0. The thermal modifications 
of hadrons are encoded in the hadron spectral functions, e.g. the presence (absence) of a resonance 
peak in the spectral function indicates the existence (dissociation) of a bound state. The hadron spectral 
function is also related to thermal dilepton and photons rates as well as transport properties of the 
medium, e.g. the electrical conductivity and heavy quark diffusion coefficient. The related hadron 
correlation functions and their susceptibilities are useful observables that allow to study the chiral 
phase transition and the U(1)a symmetry restoration. 

In the following I will give brief introductions to meson correlation and spectral functions in sec- 
tion|2], and to the chiral and U(1)a symmetries in section[3] In section|4]I will show the recent lattice 
QCD results on in-medium hadron properties. For a review on earlier results see recent reviews, e.g. 
Ref. [3]. Finally I will give a summary in section [5] 



2 Meson correlators and spectral functions 

In-medium meson properties as well as some transport coefficients are all encoded in thermal meson 
spectral functions. The meson spectral function p# (w, p) for a given quantum number channel H can 
be defined through the retarded correlator D^iio, p) 

Ph (oj, p) = D + H (<o, p) - D- H (a), p) = 2 ImD* (w, p), (1) 

where D+ (_) (w,p) = f d 4 xe iw '-^ x D + H { ~\t,x), iD+(t,x) = (J H (t, x)J H (0, 0)> and iD~ (t, x) = (J H (0,0)J H (t,x)). 
The local meson operator 7//(f, x) has the form of Jn(t, x) = tfr(t, x) ft(t, x) , with rn = 1,7s, 757^, 
for scalar {SC), pseudo scalar (PS), vector (VC) and axial vector (AV) channels, respectively. The re- 
lation of these quantum numbers to different meson states is summarized in Table [1] 
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Table 1. Meson states in different channels for light, charm and bottom quarks. 



The spectral function in the vector channel py(w, p, T) is related to the experimentally accessible 
differential cross section for thermal dilepton production, 



-p v (co,p,T), (2) 



dwd 3 p 54tt 3 to 2 (e OJ / T - 1)' 

where a em is the electromagnetic fine structure constant. The low frequency limit of the vector spec- 
tral function is related to the transport properties of the medium. In the light meson case, the spatial 
components of the vector spectral function yield the electrical conductivity 

^ = ^Umlim P '^ P ' 7 ' ) , (3) 
T 6 ^op^o u)T 

where C em is the sum of the square of the elementary charges of the quark flavor /, C em = 2/ Qf- The 
electrical conductivity quantifies the response of the QGP to electromagnetic fields and is also related 
to the emission rate of soft photons, 

dR y 3 

lim w— = —cr(T)T a em . (4) 
w^o d- 5 /? 2n l 

Furthermore, the photon emission rate of a thermal medium w^rj is also related to the light vector 
spectral function at light-like 4-momentum 

d^y_ r a em p v {uj = |p|, T) 
W dV " 47T? ePV - 1 ' ( ' 

The spatial components of the heavy quark vector spectral function are related to the heavy quark 
diffusion coefficient D 

D = -i- lim lim V H " V \ (6) 

6y°° ihOphO^J CD 

!=1 

where x°° is the quark number susceptibility that is defined through the zeroth component of the 
temporal correlator in the vector channel. The heavy quark diffusion coefficient is also related to the 
ratio of shear viscosity over entropy density tj/s H. 

On the lattice, however, spectral functions cannot be computed directly. The quantity that can be 
directly computed is the Euclidean temporal correlation function Gh(t, p). It is defined as 



? ff (r,p) = j 



d 3 xe-'> x (y H (T,x)y H (0,0)), (7) 



where Gh(t, p) is the analytic continuation of D + (t, p) fromreal to imaginary time Gh(t, p) = D + (-/t, p). 
By using the KMS relation and the above equation, one can relate the correlation function to the spec- 
tral function, 

cosh(w(T-l/2T)) 



Gh(t, p) = f ^ Ph(w, p, T) 
Jo 2 ^ 



sinh(<w/2r) 



The above equation is the basic relation for extracting spectral functions from meson correlation func- 
tions calculated in the imaginary time formalism. To extract the spectral function reliably one needs 
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large number of points in the temporal direction of the lattice. Normally the Maximum Entropy Method 
(MEM) is used to extract spectral functions from correlators [5 |. As the temporal extent, aN T , is re- 
stricted by the inverse temperature, in order to have larger number of points N T in the temporal di- 
rection one has to decrease lattice spacing a at fixed values of temperature. As a consequence the 
computing cost will increase. The Euclidean spatial correlation function Gh(z,\)±, p T ) can be com- 
puted via an integration over the Euclidean time t and two spatial directions x x = (x, y) 

G H (z, P±, Pr) = J ' dr J dx ± expHp • x) (J h (t, x)J h (0, 0)> , (9) 

where p ± = (p x ,Py), x = (x x ,t) and p = (p x ,/? T ). The spatial correlation function is related to the 
spectral function through 



C°° dp z r°° dco 2cj 

G h (z,P±,Pt) = — exp(ip z z) z r p H (uj,pj_,p z ,p T ,T). 

J-co 27T Jo 2n co z + p\ 



(10) 



The above equation is much more complicated than that of the temporal correlation function and it 
thus is not so straightforward to extract the spectral function from the spatial correlation function. 
However, one advantage of spatial correlation functions over temporal correlation functions is that the 
physical extent in the spatial direction is not restricted by the temperature and thus by going to large 
distances one can extract the exponentially decayed constant, i.e. screening mass. 

Bottomonium states are too heavy to be implemented in direct lattice QCD calculations. The lattice 
spacing required for this would need to be much smaller than the inverse of the bottom quark mass. 
A solution is provided by the Non-Relativistic QCD (NRQCD), which is an effective theory of QCD 
where physics above the scale of the heavy quark mass M is integrated out. In the framework of 
NRQCD, it is sufficient to have the lattice spacing a as coarse as M » I /a » Mu. Thus by construction 
there is no continuum limit of NRQCD on the lattice. The correlation function in NRQCD is related 
to the spectral function through 

^V(w,p,r)exp(- W r) (NRQCD). (11) 

Note that the integral kernel exp(-orr) has no temperature dependence and the correlation function 
Gh(t, p) in Eq. ( TTTb is not symmetric at r = 1/27. These are in contrast to G#(r, p) defined in Eq. |(8). 



3 Chiral and axial symmetries 

The QCD Lagrangian with Nf massless quark flavors has a large group of symmetry 

SU(Nf) L ® SU(N f ) R ® U(l) v ® U(1) A . (12) 

In the vacuum the chiral symmetry S U (N/)i ® SU (N/)r of QCD is spontaneously broken into a 
smaller SU(Nf)v vacuum symmetry. It is expected that at sufficiently high temperature the QCD 
matter will transit from a chiral asymmetric phase to a chiral restored phase. The order parameter 
for such a phase transition is the chiral condensate (tfnfr), which is nonzero below the chiral phase 
transition temperature T c and zero above T c , The chiral symmetry restoration can also be signaled by 
the degeneracy of the vector and axial vector mesons, e.g. p and a\. 

The breaking of chiral symmetry S U(3)l ® SU(3)r predicts the existence of nine massless Gold- 
stone bosons. However, there are only eight light mesons in nature, i.e. octet mesons n ± , n°, K ± , K°, 

K°, j]. The ninth "Goldstone" boson, i.e. the singlet meson 77' (958), has a much larger mass. Such a 
mass splitting is due to the violation of the axial U(1)a symmetry by the axial anomaly present at the 
quantum level [6 7 |- The effective restoration of the U(\)a symmetry may lead to a softening of the 77' 
mass which would provide interesting experimental signatures. 

The degree of U( \)a symmetry breaking may affect our understanding of QCD phase diagram. In 
the chiral limit of two flavor QCD, the order of the chiral phase transition is expected to be second 
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Fig. 1. Symmetry transformations relating pseudo scalar and scalar mesons in flavor non-singlet and singlet chan- 
nels in N f = 2 QCD @. 



order belonging to the three dimensional 0(4) universality class if the U(1)a is strongly broken near 
the phase transition. If the U{1)a symmetry breaking is negligible or U(1)a is restored near the phase 
transition then the chiral phase transition can become first order [9 |. Thus the nature of the chiral phase 
transition may depend critically on the strength of U(1)a symmetry breaking. 

The way to investigate the chiral and U(1)a symmetry on the lattice is to look at various correla- 
tion functions and susceptibilities. The susceptibilities XnAo-.n are the four volume integration of the 
corresponding correlation functions, xasc is the disconnected part of the susceptibility of the chiral 
condensate (4n//) and it equals to the disconnected part of Xa-, and^disc is the disconnected part of the 
susceptibility of (ijjy^) and it equals to disconnected part of Xn- They can also be written in terms of 
other susceptibilities as shown in Fig.Q]and thus the U(1)a and chiral symmetries are restored through 
the following degeneracies: 

U(1) A symmetry : Xn = Xs and x& +Afdisc = Xn -X5,&sa < (13) 
S U (2) L xSU (2) s symmetry : %* = Xs + *di sc and Xs = Xn ~ *5,di sc • (14) 

We see that if U(1)a symmetry is restored the values of Xn - Xs and ;f disc +^5,disc should both go to 
zero and if chiral symmetry is restored^ —Xs should equal to Xdisc and^, se — X5,disc should be zero. 
When both chiral and U(1)a symmetries are restored, ^disc snd Xs^isc should both equal to zero. 



4 Lattice QCD results 

4.1 Thermal dilepton rates, electrical conductivity and soft photon emission rate 

Recently a detailed calculation of the light vector correlation function on large isotropic quenched lat- 
tice was reported at a fixed value of temperature T - 1 .45 T c ifTOl . Calculations have been performed 
at several lattice spacings and volumes. This enables a reliable extrapolation of the vector correlation 
function to the continuum limit for a large Euclidean time interval 0.2 < tT < 0.5. From the com- 
parison of the free correlation function and the continuum extrapolated vector correlation function at 
T - 1.45 T c , it was found that the high energy part of the spectral function at T 1.45r c is quite 
close to the free field theory, thus the following ansatz for the spectral function is used to analyze the 
continuum extrapolated vector correlation function 

paid) = 2x q CBW T~7^vi +^~( l+k ) ^ tanh(w/4r) <9(« ,^U , (15) 
or + (-T/2) z 27T 

where cbw and r are fit parameters to adjust the Breit-Wigner (BW) term, and A; is a fit parameter 
accounting for corrections to the free field theory in the second term. To understand the systematic 

errors that arise from a specific fit ansatz, a function (9(wo,zl w ) = (l + e (a, cr w ~>/«H.^ j s a j so adopted in 
the second term of Eq. (TT3T >. In the limit of — > the function 0(a>o, A^) becomes a step function with 
discontinuity at wo- Fits with two extreme cases of parameter sets of ojq and A w that are constrained 
by data are shown in Fig. [2] In the left panel of Fig. [2] the thermal dilepton rates are compared to 
those calculated within the hard thermal loop (HTL) approximation [11] using a thermal quark mass 
mr/T = 1. Obviously they are in good agreement for u>/T > 2. For 1 < co/T < 2 differences 
between the HTL spectral function and lattice results is about a factor of two, which also is the intrinsic 
uncertainty in the spectral analysis. At energies u>/T < 1 the HTL results grow too rapidly, as is well 
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Fig. 2. Left: thermal dilepton rate in 2-flavor QCD. Shown are results from fits using Eq. d!5t without a cut-off 
on the continuum contribution (a>o/T = 0) and with the largest cut-off tolerable in the fit ansatz (oJq/T = 1.5). 
The HTL curve is for a thermal quark mass m T /T = 1 and the Born rate is obtained by using the free spectral 
function. Right: the spectral functions that entered the calculation of the dilepton rate i TOl . 



known. As the results shown in Fig.[2]is only for one single temperature, it is important to have results 
from lattice QCD at various temperatures around T c in the near future. 

The vector spectral function is shown in the right panel of Fig. [2] There is no resonance peak seen 
in the spectral function. Through Eq. (0 the slope of the spectral function at vanishing energy gives 
an estimate for the electrical conductivity and the zero energy limit of the thermal photon rate IfTOl 

1/3 < — - < 1, limw^ = (0.0004 - 0.0013) T 2 C at T * 1.45 T c . (16) 

C em T w-»o d s p 

Different values of the electrical conductivity from lattice QCD calculations were reported in Ref. 0121131 . 
however, these values, unlike those reported in Ref. IfTOl . are unrenormalized and not extrapolated to 
the continuum limit. The temperature dependence of the electrical conductivity has also been studied 
on the lattice recently. It was found that cr/T « 0.33^,,, and is almost independent on temperature 
in the temperature range 1.16r c < T < 2.98r c [141 . However, the results reported in Ref. [14| were 
obtained at only one value of the lattice cutoff and were not extrapolated to the continuum limit. Nev- 
ertheless these results are consistent with those in Ref. [ 10 1 . The leading order pQCD calculations 
give a much larger result, i.e. cr/T « 6 e 2 [ 15 1. Note also that the electrical conductivity obtained from 
lattice QCD is more than an order of magnitude larger than the that of a pion gas above T c ifTBI . 



4.2 Heavy quarkonia correlators and spectral functions 

The most direct approach to check whether a quarkonium state survives in the hot medium is to in- 
vestigate the temperature dependence of its spectral function. Recently a detailed analysis of the char- 
monium temporal correlation and spectral functions computed on fine and large isotropic quenched 
lattices at T — 0.73, 1.46, 2.20 and 2.93 T c was reported in Ref. fljj. Computing time can be saved 
by doing simulations on anisotropic lattice but the lattice cutoff effects are much more severe in this 
case as discussed in Ref. IfTTl . This work is mainly based on the MEM to extract the spectral func- 
tions. As the number of data points in the temporal correlation function is very important in this type 
of study, an enormous computing effort has been made to have large isotropic lattices, e.g. 128 3 x 48 at 
1.46 T c . The uncertainties on the spectral function from MEM have been carefully studied. The results 
for the charmonium spectral function in the vector channel is shown in the left panel of Fig. [3] It is 
obvious that at 1.46 T c the ground state peak seen at 0.73 T c disappears, i.e. J /iff does not survive at 
T > 1.46 T c . r] c ,XcQ and;^ c i are also found to disappear at T > 1.46 T c IfTTl . The fate of S wave sates 
(J /iff and 77 e ) suggested by the analysis in Ref. IfTTl is different from previous lattice QCD studies that 
suggested the survival of S wave sates up to 2 T c . This is mainly due to the fact that most of the previ- 
ous lattices studies were done on anisotropic lattices and had a fewer number of points in the temporal 
direction 1 17]. On the other hand, the results for P wave states (y c o and^ c i) from this work are com- 
patible with previous results, i.e. they does not survive at T > 1.46 T c . However, these results from 
Ref. [17 1 do not necessarily mean that all the charmonium states melt at the same temperature. One 
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Fig. 3. Left: Statistical uncertainties of output spectral functions in the V„ channel at temperatures T = 
0.73, 1.46, 2.20, 2.93T C . The shaded areas are statistical errors of amplitudes of output spectral functions from 
Jackknife analyses and the solid lines inside the shaded areas are mean values of spectral functions. The horizon- 
tal error bars at the first peaks of spectral functions at 0.73 T c and 1.46 T c stand for the statistical uncertainties 
of the peak location obtained from Jackknife analyses 1 17 1. Right: Pseudo scalar charmonium screening masses 
calculated with antiperiodic (filled symbols) and periodic (open symbols) boundary conditions as function of the 
temperature. Solid (dashed) lines correspond to the free theory prediction with charm quark mass m c = 1.42 GeV 
for antiperiodic (periodic) boundary conditions 1181 . 

needs further lattice QCD computations at T c <T < l.46T c to locate the dissociation temperatures of 
charmonium states. 

Spatial charmonium correlators and screening masses have also been calculated on the lattice in 
full QCD [ 1 8 1 . As shown in the right panel of Fig. [3] at T > 1 .5 T c the screening mass in the pseudo 
scalar channel differs strongly from its zero temperature value and it depends strongly on the boundary 
conditions. All this suggests the dissociation of bound states in the pseudo scalar channel at T > 1 .5 T c . 
This finding is consistent with that from the study of charmonium spectral functions extracted from 
temporal correlation functions reported in Ref . IfTTII . 



| 

3 



J __ r/r =0.42 

_ 1 — T/T =1.05 _ 

" !'* 


— r/r =1.05 




r/r =1.20 
r/r =1.40 

\ 


7/7" = 1 .40 - 

_ — r/r =1.68 - 
— ■ 


__ r/r =1.68 - 
- — r/r =1.86- 

: hl : 

|: , 1 , 1 ,V, l-^l -, 


1 


1 ' KM i- 

-- r/r =1.86 

— 77T=2.09 
\ "^(vector) 
\ Upailon 

,1,1 N-J-T-L, 



to [GeV] 

Fig. 4. Bottomonium spectral functions in the vector channel from Non-Relativisitic QCD calculations 1 19 1. 

For the fate of bottomonia at finite temperature, a lattice calculation in the framework of NRQCD 
has been performed recently 1 19]. Lattice QCD simulations were performed with two light quarks on 
highly anisotropic lattice using NRQCD for bottom quarks. The spectral functions shown in Fig. |4] 
are extracted from correlations functions using the MEM. It can clearly be seen from that figure that 
the ground state in the T channel survives up to 2.09 T c . Results for the rjb channel are similar. The 
calculation has also been extrapolated exploringly to nonzero baryon density [20|. However, one has 
to note that there is no continuum limit for the lattice version of NRQCD and the calculations in 
Ref. [19 20 1 were performed on highly anisotropic lattices that can distort the physics. 

The properties of heavy quarkonia in the vacuum have been studied based on a non relativistic 
Schrodinger equation using either modeled or lattice QCD computed heavy quark potentials. However, 
the definition of the heavy quark potential (free energy or internal energy ) is ambiguous at nonzero 
temperature. Recently it has been proposed by authors in Ref. [21 ] that in the nonrelativistic limit the 
leading order potential in an effective Schrodinger equation can be appropriately defined by assuming 
that the spectral density of the Wilson loop has peaks according to Breit-Wigner distributions. In this 
framework, the peak location of the spectral density corresponds to the real part of the heavy quark 



Hadronic resonance production in heavy ion and elementary collisions 



potential while the peak width corresponds to the imaginary part of the heavy quark potential. Using 
the MEM it is feasible to extract the peak location, however, it is very hard to extract the width of the 
peak reliably. 

The charm diffusion coefficient obtained from the slope of the vector spectral function is almost 
XjnT at T e (1.46 T e ,2.93 T c ) and compatible with zero below T c ("171 . The mean value of the 
heavy quark diffusion coefficient obtained from the computation of heavy quark effective theory on the 
lattice is somewhat larger B23I241 . i.e. 2/tzT but compatible with the charm diffusion coefficient from 
Ref. [ 17 1 within errors. These values from lattice QCD are smaller than that from pQCD calculations in 
the investigated temperature window 11251 and lies in the range used in phenomenological models 1 26 1 . 



4.3 Chiral phase transition and U(l) A symmetry restoration 

Studies on this topic have also been carried out recently in Ref. ET\ . In this study the improved (p4) 
staggered action was used and the effective restoration of the U(1)a symmetry occurs at a temperature 
larger than the QCD transition temperature indicated from the degeneracy of the screening masses 
of pseudo scalar and scalar correlators. A better discretized version of staggered fermions has also 
been implemented in Ref. [28 1 to investigate the fate of the axial symmetry at finite temperature and 
their results suggests that U{1)a symmetry still remains broken at about 1.1 T c . Due to subtle issues 
concerning the realization of chiral and axial symmetries in the staggered discretized fermion action, 
better discretization schemes, e.g. domain wall fermions and overlap fermions, which preserve the full 
chiral symmetry of continuum QCD and reproduce the correct axial anomaly even at nonzero lattice 
spacing, have also been implemented to study this topic. 



150 160 170 
T (MeV) 




Fig. 5. Left: susceptibilities as function of temperature, Right: the time histories of the topological charge Q top 
and;^ —Xs< Figures are taken from Ref. (8). 



Using domain wall fermions the HotQCD collaboration recently investigated the restoration of the 
spontaneously broken chiral symmetry and the effective restoration of the anomalously broken U(1)a 
symmetry at finite temperature on 16 3 x 8 lattices with m„ « 200 MeV [8]. One of the main results 
is shown in Fig. [5] The three quantities, x disc, A's.disc and;^ ~Xs> agree within errors for temperature 
T > 170 MeV. This suggests the restoration of the S U{2)i x S U(2)r symmetry. The U(1)a symmetry, 
on the other hand, remains broken up to temperature T * 200 MeV as indicated from the nonzero 
value of^ - X&- If F70 MeV is taken as the pseudo critical transition temperature, this indicates that 
U(1)a still remains broken at T ~ l.2T pc . As shown in the right panel of Fig. [5] the nonzero value 
°f Xn ~ Xs is related to the nonzero topological charge Q lop , which will be suppressed in the large 
volume limit. To draw reliable conclusion on U(1)a symmetry restoration at high temperature, one 
needs to carry out simulations with larger volumes and perform detailed analysis on the Dirac opera- 
tor eigenvalue spectrum. The JLQCD collaboration using overlap fermions performed simulations at 
higher temperatures mainly on 16 3 x 8 lattices with m„ > 290 MeV and Q top fixed to be zero [29|. 
They found that U(1)a symmetry is restored at 209 MeV by looking at the degeneracies of correlators 
and the gap around origin in the Dirac operator eigenvalue spectrum. However, the volume effects on 
topological charge fixing might be large and need to be checked. 

From the studies discussed above, the effective U(1)a symmetry restoration seems to occur after 
the chiral symmetry restoration and U(1)a restoration temperature certainly needs further studies. Of 
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course, it would be very interesting to have good signals to probe the possible restoration of U(1)a 
symmetry in the heavy ion experiments 113013 11321 . 



5 Summary 

The lattice QCD calculations presented here show that there has been significant progress in our un- 
derstanding of strongly interacting matter through the in-medium hadron properties in recent years. 
However, more needs to be done on the lattice to help to understand the heavy ion experiment data. 
The U(1)a symmetry restoration that is being intensively studied on the lattice may be observed in the 
heavy ion experiments from some interesting signals, e.g. a softening of the r( mass. 
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